Multiple flux jumps and irreversible behavior of thin Al superconducting rings 
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An experimental and theoretical investigation was made of flux jumps and irreversible magnetiza- 
tion curves of mesoscopic Al superconducting rings. In the small magnetic field region the change of 
vorticity with magnetic field can be larger than unity. This behavior is connected with the existence 
of several metastable states of different vorticity. The intentional introduction of a defect in the 
ring has a large effect on the size of the flux jumps. Calculations based on the time-dependent 
Ginzburg-Landau model allows us to explain the experimental results semi-quantitatively. 
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I. INTRODUCTION 

Recently, Pedersen et aL0 observed jumps in the mag- 
netization of superconducting rings which corresponds to 
changes of the vorticity larger than unity. This is in 
contrast to the behavior of superconducting disks where 
only changes in the value of the vorticity of unit size 
were observeda In some respect the observed behavior 
in rings is similar to vortex avalanches which werjEppb- 
served in superconductors with strong bulk pinningj'u or 
to jumps in the magnetization when several vortices (in 
the form of a chain) enter in a supercondujcting film of 
width comparable to the coherence lengthOEl. The occur- 
rence of such jumps in a defect free superconducting ring 
originates from the fact that several metastable states 
with different vorticity L are possible for a given mag- 
netic field. However the existence of such multiple stable 
states is not a sufficient condition to explain changes in 
the vorticity larger than unity (e.g. they also exist in the 
case of superconducting disks). An additional important 
requirement is to find the stability condition and finally 
the state to which the system relaxes to. This requires 
the study of the transition process from one state to an- 
other, i.e. it requires the analysing the time-dependent 
process. 

The stability condition was studied numerically in 
Ref.Q for the case of a hollow cylinder, and in a number 
of works (see for example Ref.Q'tl and references therein) 
for superconducting disks and rings by using the static 
Ginzburg-Landau (GL) equations. Unfortunately no an- 
alytical results were presented due to the rather general 
character of the studied systems in the above works. 

Recently, we studied the transition procesaliJ using 
the time-dependent Ginzburg-Landau equations. It 
was shown that transitions between different metastable 
states in a mesoscopic superconducting ring are governed 
by the ratio between the time relaxation of the phase of 
the order parameter (which is inversely proportional to 



the Josephson frequency) and the time relaxation of the 
absolute value of the order parameter r|^| . We found that 
if the ratio t\^\/t^ is sufficiently large the system will al- 
ways transit from a metastable state to the ground state. 
This leads to an avalanche-type variation of L when the 
vorticity of the metastable state differs appreciably from 
the vorticity of the ground state. In contrast to the case 
of a superconducting film, in a ring the 'vortex' entry 
occurs through a single point and the vorticity increases 
one by one during the transition. In low-temperature 
superconductors like In, Al, Sn the ratio t\^\It^ is very 
large for temperatures far below the critical temperature 
Tc and hence, if such systems are driven far out of equi- 
librium they will always relax to the ground state. 

In this work we investigate the conditions under which 
a state with a given vorticity becomes unstable in a finite 
width ring and we find how the supeconducting order pa- 
rameter in the ring changes with increasing applied mag- 
netic field. We are able to find an analytical expression 
for the dependence of the order parameter on applied 
magnetic field, and hence for the upper critical field at 
which superconductivity vanishes in such a sample. We 
provide a direct comparison between the theoretical and 
experimental results on aluminium rings. Our theoreti- 
cal calculations are based on a numerical solution of the 
time-dependent Ginzburg-Landau equations. 

The paper is organized as following. In section II 
the theoretical formalism is presented and the two- 
dimensional time-dependent GL equations are solved. In 
Section III the experimental results are presented and 
compared with our theory. In Section IV we present our 
conclusions and our main results. 



II. THEORY 

We consider sufficient narrow rings such that we can 
neglect the screening effects. This is allowed when the 
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width of the ring w is less than maxX, X"^ /d, where A is 
the London penetration length and d is the thickness of 
the ring. In order to study the response of such a ring 
on the applied magnetic field we use the time-dependent 
Ginzburg-Landau equations 



. (V-iA)V+(l-|^nV', (la) 
Aip = div(Im(V'*(V-iA)?A)), (lb) 



where all the physical quantities (order parameter -0 = 
IV'le*'^, vector potential A and electrostatical potential 
ip) are measured in dimensionless units: the vector po- 
tential A is scaled in units <i>o/(27r^) (where $0 is the 
quantum of magnetic flux), and the coordinates are in 
units of the coherence length £,{T). In these units the 
magnetic field is scaled by Hc2 and the current den- 
sity, j, by jo — c$o/87r^A^^. Time is scaled in units of 
the Ginzburg-Landau relaxation time tqe = 47r(T„A^/c^, 
the electrostatic potential, ip, is in units of c$o/87r^^A(j„ 
(cr„ is the normal-state conductivity). Here the time- 
derivative is explicitly included which allows us to deter- 
mine the moment at which the state with given vorticity 
L becomes unstable. It is essential to include the electro- 
static potential (which is responsible for the appearance 
of the Josephson time or frequency) in order to take into 
account the multi-vortex^Jiajps. In some previous stud- 
ies (see for example Refs.EJ't5) ip = was assumed and as 
a consequence only transitions with unit vorticity jumps, 
i.e. AL = 1, are possible in the ringtj. The coefficient 
M = 48 was chosen such that after the transition tbe sys- 
tem is in the thermodynamically equilibrium stateli3. We 
assume that the width (w) of the ring is less than two 
coherence length f , because: i) all experimental results 
presented here were performed for such type of samples; 
and ii) only in this case it is possible to obtain simple 
analytical expressions. For instance, the dependence of 
the order parameter on the applied magnetic field and 
the upper critical field Hmax- 

For w < 2^ the order parameter is practically indepen- 
dent of the radial coordinate. This is demonstrated in 
Fig. 1 where the dependence of the order parameter in 
the middle of the ring is compared with its value at the 
inner and outer boundary of the ring, i.e. r = RzLw/2 
{R is the mean radii of the ring), for two different rings. 
Notice that these two numerical examples corresponds al- 
ready to relative thick mesoscopic rings, i.e. R/w ^ 1 — 2. 
For the field Hmax we are able to fit our numerical results 
to the expression 



H„ 



3.67- 



(2) 



For rings with w < 2^ and w/R < 1 this analytical ex- 
pression is within 2% of the numerical results. It is inter- 
esting to note that Hmax does not depend on the radii 
of the ring. But the value of the vorticity of the sys- 
tem depends on R. For example, for R = 5.5^(16.5^) 
L = 55(501) for w = ^ at = Hmax- 
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FIG. 1: Dependence of the absolute value of the order pa- 
rameter 1 7/) I on the applied magnetic field for two different 
rings in the ground state. Dashed curve corresponds to 
\tj}\{R — ■w/2,H), solid curve to \^\{R,H) and dotted curve 
to \rp\{R + w/2). 



Note that Eq. (2) has the same dependence on the su- 
perconducting parameters as in the case of a thin plate 
with thiickness d < \/5A placed in a parallel magnetic 
fieldlljo. Even the numerical coefficient is quite close, 
i.e. for a thin plate it is equal to 2\/3 ~ 3.46. Further- 
more, we found that the transition to the normal state 
of our rings at the critical field Hmax is of second order 
as is also the case for a thin plate. A possible reason for 
this close similarity is that for a thin plate with thick- 
ness d < \/5X the screeniag|-cffects are also very small. 
In the calculations of Refs.E3o an averaged value for the 
order parameter was used independent of the coordinate. 
Note that this is similar to our \tp\ which is practically 
independent on the radial coordinate (see Fig. 1). 

The absolute value of the order parameter (in the mid- 
dle of the ring) is, too a high accuracy, given by the ex- 
pression 



\lP\^ = 1-{H/Hmax?-P{L,H) 



(3) 
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with p{L, H) = L/R — HR/2, where the vorticity L de- 
pends on the history of the S5pteHi. This result is similar 
to the one obtained in Refs.tjEJ with the exception of 
the last term in Eq. (3) which appears due to the closed 
geometry of the ring and hence leads to a nonzero L. 
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FIG. 2: Magnetic field dependence of the magnetization (a), 
the order parameter (b) and the gauge invariant momentum 
(c) in the middle of the ring. Dotted curve in Fig. 2(b), is the 
expression VI - {H/Hma.y. Dotted curve in Fig. 2(c) is the 
expression ^l-{iH-Ho)/H„,a.)yV3, where Ho ~ 13G - 
is the displacement of the maximum of M{H) from the H = 
line. 

All the above results were obtained for a ring which is 
in the ground state at any value of the magnetic field. 
However such a system can exhibit several metastable 
states at a given magnetic field, and consequently this 
may lead to hysteretic behavior when one sweeps the 
magnetic field up and down. Furthermore, with changing 
field the vorticity may jump with AL > I. An impor- 
tant question which arises is the condition of stability of 
the state with given vorticity. This-cmestion was stud- 
ied earlier for one-dimensional ringsE3'E3, i.e. rings with 
zero width. It turns out that the system transits to a 
state with another vorticity when the value of the gauge- 



invariant momentum p 
value 



V(f) — A reaches the critical 



(4) 



At this condition it is easy to find the value of the field 
for the first 'vortex' entry 



V^R 



1 

2R^' 



(5) 



We will now generalize the results of RefsEjO to the 
case of finite width rings with w < 2^. First we will 
neglect the dependence of ip on the radial coordinate in 
which case the GL equations reduce to one-dimensional 
expressions. But in order to include the suppression of 
the order parameter by an external field for a finite width 
ring we add the term —{H/Hmax)^^ to the RHS of Eq. 
(fa), where Hmax is given by Eq. (2). Using the stabil- 
ity analysis of the linearized Ginzburg-Landau equations 
near a specific metastable state as presented in RefJl3 we 
obtain the modified critical momentum 



(6) 



Note that now Pc decreases with increasing magnetic 
field. This automatically leads to a decreasing value of 
the jump in thje vorticity AL at high magnetic field, be- 
cause m Ref.B it was shown that 




ALraax = Nint(]3c-R), 



(7) 



where Nint(a;) returns the nearest integer to the argu- 
ment. 

In order to check the validity of Eq. (6) we performed 
a numerical simulation of the two-dimensional Ginzburg- 
Landau equations, Eqs. (fa,b), for a ring with R — 5.5^ 
and w = 1.5^ (for these parameters the theoretical find- 
ings fit the experimental results - see section below). In 
Fig. 2 the magnetization, the order parameter and the 
gauge-invariant momentum p are shown as function of 
the applied magnetic field. The magnetic field was cy- 
cled up and down from H < —Hmax to H > Hmax- The 
condition (6) leads to an hysteresis of M{H) and to a 
changing value of the jump in the vorticity in accordance 
with the change in pc- The main difference between our 
theoretical prediction (6) and the results of our numerical 
calculations appears at fields close to Hmax- Apparently 
it is connected with the fact that for the considered ring 
the distribution of the order parameter along the width 
of the ring is appreciably nonuniform at ~ Hmax and 
as a consequence the one-dimensional model breaks down 
(see Fig. I). 

Finally, we also considered the same ring with a defect. 
The effect of the defect was modelled by introducing in 
the RHS of Eq. (fa) the term —p{s)'4> [s is the arc- 
coordinate) where p{s) = —1 inside the defect region 
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FIG. 3: Magnetic field dependence of the magnetization (a), 
the order parameter (b) and the gauge invariant momentum 
(c) (in the middle of the ring) of a ring containing one defect. 
Dotted line in Fig. 3(b) is the expression \/l — (H/Hmax)^ ■ 



with size ^ and p{s) = outside. This leads to the results 
shown in Fig. 3 for M(i7), |-0|(i/) and p{H). Due to the 
presence of the defect, Pc differs from Eq. (6) already at 
low magnetic field {pc{H = 0) ~ 0.33 at given "strength" 
of the defect) and as a result only jumps with AL = 1 are 
possible in such a ring. In this case the Pc and also 
depend on the applied magnetic field with practically the 
same functional dependence on H as Eq. (6). 



III. COMPARISON WITH EXPERIMENT 

The measurements were performed on individual 
Al superconducting^pdngs by using ballistic Hall 
micromagnetometrjOO. The techniques employs small 
Hall probes microfabricated from a high-mobility two- 
dimensional electron gas (2DEG). The rings - having 
radii R ~ 1/im and width w ranging from 0.1 to O.Sfim 
- were placed directly on top of the microfabricated Hall 
crosses, which had approximately the same width b of 
about 2fj,m (see micrograph in Fig. 4 for a ring with 
an artificial defect). These experimental structures were 



prepared by multi-stage electron-beam lithography with 
the accuracy of alignment between the stages better than 
lOOnm. The rings studied in this work were thermally 
evaporated and exhibited a superconducting transition at 
about 1.25K. The superconducting coherence length was 
f (T = 0) ~ 0.18/im. The latter was calculated from the 
electron mean free path / ~ 2bnm of macroscopic Al films 
evaporated simultaneously with the Al rings. The Hall 
response, Rxy^ of a ballistic cross is given by the amount 
of magnetic flux J-P4s through the central square area 
6 X 6 of the crosaijlla. For simplicity, one can view the 
ballistic magnetometer as an analogue of a micro-SQUID, 
which would have a square pick-up loop of size b and 
superconducting rings placed in its centre. We present 
our experimental data in terms of the area magnetiza- 
tion M =< B > —H which is the difference between the 
applied field H and the measured field < B >^ Rxy 
Previously, we have studied individual superconducting 
and ferromagnetic disks apd found excellent agreement 
with the above formulaEZlEj. For further, ^ftflilf about 
the technique, we refer to our earlier worklijOo. 




FIG. 4: A micrograph of the superconducting ring placed on 
top of a Hall bar. An artificial defect (narrowing of the ring 
cross section) is intentionally made by electron beam lithog- 
raphy. 

Rings with and without an artificial defect were stud- 
ied. Let us consider first the ring without artificial defect. 
In Fig. 5(a) the full magnetization loop of such a ring 
with parameters R — 1.0±0.1/im and w — 0.25±0.05/ir7i 
is shown. In Fig. 6 (solid curve) the low field part 
of the virgin curve is presented. From the virgin trace 
M{H) we can find the magnetic field for the first vor- 
tex entry, H^n, and hence we estimate ^ ~ 0. 19 fim at 
the given temperature (T ~ OAK) using Eq. (5) (this 
value of ^ is in agreement with the above experimental 
value 1^(0) ~ 0.18/im obtained from the mean free path). 
Furthermore, we know from Fig. 5(a) that the vorticity 
changes with AL = 3 for ~ 0. This agrees with the 
fact that the radii of the ring is larger than 4.6^ (see Eq. 
(7)). Another important information which may be ex- 
tracted from the virgin curve is that at the first vortex 
entry the magnetization drops considerably but it does 
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not change sign. If we recall that at every vortex entry 
p decreases on 1/R (and hence the current density j and 
M ~ / [j X Y\dV also changes proportionally) we can con- 
clude that the radii of our ring should be in the range 
5.5^ < R < 6.5^. This agrees with the experimental 
value R/^~ 5.3 ±0.5. 
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FIG. 5: Magnetic field dependence of the magnetization of the 
ring without (a) and with (b) an artificial defect at T ~ OAK. 
Parameters of the rings (width and radii) are the same within 
experimental accuracy. 

If we take the above value for ^ and w ~ 1.5^ we ob- 
tain the maximum field of Hmax — 223 G. This value 
is slightly smaller than the value obtained from Figs. 2 
and 3 Hmax — 2400 which we attribute to the large 
coordinate step which we used in our numerical calcu- 
lations of Eq. l(a,b). The value is also larger than the 
experimental value Hmax — 185 G. This disagreement 
between theory and experiment is most likely connected 
to the semi-quantitative applicability of the Ginzburg- 
Landau equations in the considered temperature range. 
The range of applicability of the Ginzburg-Landau equa- 
tions (even stationary ones) for this specific supercon- 
ductor is very narrow. Nevertheless based on previous 
comparison between experinteote and theory for meso- 
scopic superconducting disksE3Eil it was found that the 
GL equations provided a rather good description of the 



superconducting state even deep inside the (iJ, T) phase 
diagram. 

Figs. 2(a) and 5(a) are qualitatively very similar. For 
example our theory describes: i) the hysteresis; ii) the 
non-unity of the vorticity jumps, i.e. AL = 3 in the 
low magnetic field region, AL = 2 in the intermediate 
H- region, and ALT = 1 near Hmax- Theoretically (ex- 
perimentally) we found 6(5), 13(21), 22(18) jumps with 
respectively AL = 3, 2, 1; and iii) the non symmetric 
magnetization near ±Hmax for magnetic field sweep up 
and down. 

In the ring with approximately the same mean radii 
and width but containing an intentionally introduced ar- 
tificial defect, jumps with AL — 1 arc mostly observed 
(see Fig. 5(b)). The reason is that an artificial defect 
considerably decreases the critical value Pc (and hence 
the field Hen - see dotted curve in Fig. 6). From Figs. 
2(c), 3(c) it is clear that the maximum value p^/ — 0.54 
for a ring without defect and pj? ~ 0.35 for a ring with 
a defect. The ratio p^/p^^ — 0.65 is close to the ratio of 
the field of first vortex entry H^^^/Hlf^ ~ 0.67 obtained 
from experiment (see Fig. 6). From Fig. 2(c) it is easy 
to see that for a ring without a defect at p ~ 0.35 there 
are only jumps with AL = 1. But if we slightly increase 
p then jumps with AL ~ 2 can appear in the system. 
So we can conclude that p — 0.35 is close to the border 
value which separates regimes with jumps in vorticity of 
AL = 1 and AL = 2. From our experimental data it fol- 
lows that the maximum value of Pc is very close to this 
border. Thermal fluctuations may influence the value of 
AL, in particular for a,pc value close to this border value. 
This is probably the reason why in the experiment (Fig. 
5(b)) occasional jumps with AL = 2 are observed which 
are absent in our simulation (Fig. 3(a)). 
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FIG. 6: Magnetic field dependence of the virgin magnetization 
of a ring without (solid curve) and with (dotted curve) an 
artificial defect. The dotted curve is shifted for clarity by 0.6 
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IV. CONCLUSION 

We studied multiple flux jumps and irreversible behav- 
ior of the magnetization of thin mesoscopic Al supercon- 
ducting rings. We have shown experimentally and the- 
oretically that at low magnetic fields and for rings with 
sufflciently large radii the vorticity may change by val- 
ues larger than unity. With increasing magnetic field the 
order parameter gradually decreases and thus leads to a 
decrease of the size of the jumps in the vorticity. For 
rings with width less than 2^ analytical expressions were 
obtained for the dependence of the order parameter on 
the applied magnetic field. We have found that a state 
with a given vorticity becomes unstable when the value 
of the gauge-invariant momentum reaches a critical value 
Pc which decreases with increasing magnetic field. This is 
responsible for the fact that AL decreases with increas- 



ing H. The introduction of an artificial defect in the ring 
leads to a decrease of Pc in comparison to the case of a 
ring without a defect and also results in a decrease of 
AL. 
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